Reconstructing supernova-neutrino spectra using low-energy beta-beams 
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Only weakly interacting, neutrinos are the principal messengers reaching us from the center of 
a supernova. Terrestrial neutrino telescopes, such as SNO and SuperKamiokande, can provide 
precious information about the processes in the core of a collapsing and exploding star. But the 
information about the supernova that a neutrino detector can supply, is restricted by the fact that 
little experimental data on the neutrino-nucleus cross sections exists and by the uncertainties in 
theoretical calculations. In this letter, we propose a novel procedure that determines the response of 
a target nucleus in a supernova-neutrino detector directly, by using low-energy beta-beams. We show 
that fitting 'synthetic' spectra, constructed by taking linear combinations of beta-beam spectra, 
to the original supernova-neutrino spectra reproduces the folded differential cross sections very 
accurately. Comparing the response in a terrestrial detector to these synthetic responses provides a 
direct way to determine the main parameters of the supernova-neutrino energy-distribution. 

PACS numbers: 25.30.Pt,26.50.+x,97.60.Bw,13.15.-|-g 



Understanding the mechanism of a Type-II core- 
collapse supernova is a longstanding problem. While the 
general nature of these events is understood to be the ex- 
plosion that results from the collapse of a massive star, 
a completely self-consistent model of the explosion does 
not yet exist. Weak physics and neutrino-interactions 
are widely regarded as central players in the explosion 
mechanism. The core is so hot and dense that even the 
neutrinos are trapped. These neutrinos scatter out of the 
core on a timescale of about ten seconds, taking away 
with them the vast majority of the gravitational binding 
energy released in the collapse. The handful of neutrinos 
that were detected from Supernova 1987A confirmed this 
basic picture of the nature of core-collapse supernovae. 

The development of neutrino facilities and telescopes 
during the last decades ensures that a similar event today 
would provide much more detailed information. When 
the next galactic supernova occurs, neutrino detectors 
will be on the watch, using different targets and detec- 
tion mechanisms to catch neutrinos of various energies 
and flavors, thus gathering information about supernova 
(SN) neutrinos and about the supernova processes that 
generated and released them. Since the neutrinos scatter 
last at the neutrino sphere within the stellar core, their 
energy spectrum encodes precious information about the 
center of these astrophysical phenomena. 

The extraction of the information SN neutrinos are car- 
rying, however faces several intertwined problems. Un- 
certainties in the description of neutrino scattering with 
nuclear matter, general relativistic effects, and the ef- 
fect of rotation are carried over into the related produc- 
tion and release of neutrinos, and induce uncertainties in 
the energy spectra and flavor distributions of the emit- 
ted particles [ij, l^i b| - Neutrino oscillations will mix the 
original spectra on their way out of the exploding star, 
thereby veiling the information they can provide jj,|l 



Beyond doubt, unraveling the information SN neutrinos 
deposit in a detector involves the need for an optimiza- 
tion of the interpretation of the signal. 

Neutrino-nucleus scattering is an important detection 
mechanism. The cross sections are relatively large and 
reactions occur in the charged (CC) as well as neutral- 
current (NC) channels. The main nuclei of interest are 
^^O in SuperKamiokande jg, the deuteron for SNO [y], 
s^Fe in LVD Q, ^Spe and ^ospb for the proposed LAND 

iB and OMNIS [ij detectors and ^^C for KaniLAND 

3, BooNe 0, and LENA |li|. 

The extent to which we can determine the physics 
of the supernova from what we see in the detector de- 
pends on our understanding of the neutrino- nucleus cross 
sections. However, the available experimental neutrino- 
nucleus scattering data in the relevant energy range is 
limited. Neutrino cross sections are small and monochro- 
matic neutrino beams are not available. The shape of the 
neutrino energy-distribution experimentally in use thus 
far [l^ |la| , while in the same energy range as that of the 
SN neutrinos, is not the same as an SNi/ spectrum. More- 
over, theoretical estimates of cross sections and response 
details at these low energies are largely subject to nuclear 
structure uncertainties and model dependencies. Fig. ^ 
shows that the uncertainties in ^^^Pb{v, v') cross section 
calculations of '^^^Pb{v, v') induce large deviations in the 
corresponding folded SNz/ responses. The overlap of the 
shaded regions in the figure make clear how uncertainties 
in cross sections calculations will impede the reconstruc- 
tion of the original spectrum's parameters from an SNi/ 
signal. It is these problems that we want to address using 
the strength and flexibility of the low-energy beta-beam 
concept |2^ . Low-energy beta-beam facilities would pro- 
vide intense and well-characterized neutrino beams with 
energies typically in the 10 to 100 MeV range [2llEi|. 
In this letter, we show that taking appropriate combi- 
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FIG. 1: Uncertainty on theoretical cross sections for NC 
neutrino-scattering ofT ^"^^Pb, folded with different PL spec- 
tra, characterized by average energy (e^) and width parameter 
a ^3- The cross sections were taken from Refs. P. fia. Il9|l . 
The full line represents the average of these predictions, the 
borders of the shaded regions are determined by the largest 
and smallest values of the results 0, Ha, [l3 • 



nations of beta-beam neutrino cross-sections allows one 
to reconstruct differential SNz^ cross-sections, and that 
this technique provides an efficient way to reconstruct the 
SNi/ spectrum from the signal it induced in a detector. 
This solves two problems : first, we avoid uncertainties 
and model dependencies rising in theoretical cross-section 
predictions. Second, we address the complications raised 
by the lack of monochromatic neutrino beams. 

Since (anti)neutrinos of different energies and flavors 
decouple at different points in the core, the neutrino 
spectrum nsNi^v) is not a purely thermal one. Com- 
mon parameterizations are given by a Fermi-Dirac (FD) 
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jTT i,^E^/T-r, or a power- law (PL) spectrum 
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tions yield similar energy distributions, characterized by 
the average energy or temperature T of the neutrinos 
and the width of the spectrum. The average energies 
{Eu) vary roughly between 12 and 18 MeV for v^ and 
i^e's, and between 16 and 22 MeV for v^ and Vr- The 
width of the spectrum characterized by the parameters rj 
or a influences the weight of its high-energy tail |l7l | . 

Similar to the shape of SNi/ spectra, beta-beam spec- 
tra are characterized by long tails, while their average 
energy and precise shape depend on the boost-factor 7 of 
the primary beam and the geometry of the experimental 
setup [20, I21I I23I ■ The folded cross sections for neutrino- 
scattering processes show a very smooth behavior as a 
function of the 7- factor of the beta-beam spectrum |21| . 
Hence, a limited number of measurements can be used 
to infer the full 7-dependence of the folded cross section. 
Here, we suggest representing the SNi^ spectra n.siy as 
linear combinations of beta-beam spectra n^le^) J23l |. by 
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The energy transfer to the target nucleus and the ex- 
cited state uj it is left in, determine the decay mode of 
the nucleus. Therefore, the reaction products that will 
be observed in the detector are fully determined by the 
differential folded response. In order to keep the discus- 
sion as general as possible, we do not restrict our study 
to e.g. electron emission or nucleon knockout, but we 
will use full differential cross sections to illustrate the 
proposed procedure. 

The detector response to a beta-beam neutrino spec- 
trum is then given by J ds^a{e^,Lu)nr'^{e^). An actual 
supernova will produce a neutrino response determined 

by crf°g^ai('^) = liT de^aie,,, uj)nsrgnai{£,.)- Adopting 
the philosophy of Eq. ^ this quantity can be represented 
as a linear combination of beta-beam responses : 
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Once a measured SN response has been obtained, one 
can find the expansion parameters a^'^ , 7j by fitting the 
neutrino signal to the combinations of Eq. El In the re- 
mainder of the paper, we first show that only a few beta- 
beam spectra in the linear combination Q are sufficient 
for an accurate reconstruction of the SNz^ signal, then we 
illustrate that using this technique it is possible to deter- 
mine the neutrino spectrum from the supernova without 
knowledge of cr(e, w). 

The reconstruction is realized by minimizing the ex- 
pression J de \nfit(£v) — nsN(sv)\, varying the expansion 
coefficients a^' and boost factors 7^. In Fig. |2]we show 
the corresponding differential folded cross sections for an 
^^O target and a typical SNi/ spectrum. It is clear that 
the total strength as well as the position and width of the 
resonances are reproduced in a very convincing way, for 
a very small number N of beta-beam spectra in the fit. 
The method is equally successful in reproducing the sep- 
arate multipole contributions. This provides evidence for 
the usefulness of the proposed method in predicting the 
yield in an SNi/ detector, without going through the in- 
termediate step of using a nuclear-structure calculation. 
The deuteron has a lower threshold than 'massive' nuclei, 
and cross sections peak at substantially lower energies. 
This has as consequence that constructed spectra must 
contain a beta-beam component at fairly low 7. Fig. |21 
shows that the fit works well in reproducing the position 
and strength of the resonances. For spectra with smaller 
average energies the agreement is improving considerably 
when including lower gammas in the fit. Fig. 01 illus- 
trates that using the proposed technique, the agreement 
between original and reconstructed response is sufficient 



^ 0.1 



'0.01 




FIG. 2: Comparison between differential cross section for NC 
scattering on O p^ . folded with an FD distribution with 
average energy 18 MeV and width, characterized by the de- 
generacy parameter, of 77=1.1 [13, and synthetic spectra with 
N=3 (dotted) and N=5 components (dashed). The choice of 
7 was restricted to integer values between 5 and 15. The ma- 
jor multipole contributions (long-dashed) are compared with 
results of the A'^ = 5 fit (dashed-dotted). 
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FIG. 3: Comparison between differential cross sections ',24| for 
the CC reaction Ve + d—^p+p + e', folded with a PL SN;/ 
spectrum (full line) and synthetic spectra with 5 components, 
using 7-factors ranging from 5 to 15 (dotted) and including 
lower energy beta-beam spectra with 7=8 to 15 (dashed). 



to discriminate between spectra with only small separa- 
tions in average energy. The right panel indicates that 
the method is able to distinguish between cross-section 
differences induced by more subtle variations in the spec- 
tra, such as slightly different widths. 

We now turn to our technique for extracting informa- 
tion about the SN neutrinos by combining the response in 
the SN detector and information from low-energy beta- 
beams. That this is not a trivial task, is demonstrated 
in Fig. El showing that different SNi/ spectra can be ex- 
tremely similar, e.g. when having higher average energies 
but a reduced width. Obviously an efficient tool to dis- 
tinguish between these spectra is needed. 



Given the SNv signal and a set of beta-beam response 
data, fitting the SNi/ response in the detector to a lin- 
ear combinations of beta-beam responses by minimizing 
will yield the N expansion parameters 



fold 

^fit 



a 



fold 
signal I 



a^' and 7-values 7* that offer best agreement with the 
supernova data. Eq.^then immediately reveals the orig- 
inal SNi/ energy-distribution. The agreement in the pre- 
dicted responses obtained for remarkably small values of 
N guarantees the efficiency of this formalism. The ac- 
curacy of the technique is restricted mainly by detector 
limitations and is basically free of theoretical ambiguities. 

As a full Monte-Carlo analysis of supernova and beta- 
beam responses is as yet untimely, we illustrate that the 
proposed procedure can successfully be inverted and that 
it is stable against experimental uncertainties, using a 
schematic model. We create synthetic data represented 
by a cross section, folded with a synthetic neutrino energy 
distribution. The original SNiy spectrum nsN that offers 
the best fit to these synthetic data is then determined 
by selecting the spectrum that minimizes the difference 
between a member of the cr^^ set and the synthetic re- 
sponse. The results of this process can be seen in Fig. El 
In view of the problems suggested by Fig. \^ the agree- 
ment is very good, especially at higher energies. Experi- 
mental uncertainties in the beta-beam measurement and 
detector signal will give rise to noise on the expansion 
parameters that are determined from fitting the SNz^ sig- 
nal to a combination of beta-beam data. To estimate the 
effect of this on the reconstruction of the SNi/ spectrum, 
we repeated the procedure for values of the coefficients Oi 
which vary up to 10%. The mild scatter this produces, 
shows that the method is stable against uncertainties at 
this level. 

In practice, the signal in the detector will be the su- 
perposition of the Ve, ^e and the heavy-flavor spectra. 
Oscillations will mix these original spectra. Though that 
leaves the NC signal unchanged, the flavor shifts between 
heavy and electron (anti)neutrinos will most likely induce 
linear combinations of the undistorted CC response. In 
terms of the proposed procedure this means that one has 
to deal with an 'extra' linear combination superimposed 
on the one of Eq. ^ Comparing the actual combination 
observed in NC and CC channels will provide information 
about the way the neutrinos were mixed by oscillations. 

Summarizing, we studied the potential use of a low- 
energy beta-beam facility for the understanding of su- 
pernova physics. We showed that using appropriate lin- 
ear combinations of beta-beam spectra allows for a very 
accurate reproduction of supernova-neutrino responses. 
The proposed procedure allows one to make the most 
efficient use of the neutrino signal to filter information 
about the supernova from the response in a terrestrial 
detector, thus showing that beta-beams are able to pro- 
vide information about the spectra supernova-neutrinos 
were released with at the explosion site, in a direct way. 
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FIG. 4: Comparison between cumulative cross sections for NC scattering on ^"^Pb, folded with a PL SNz/ spectrum (full line) 
and synthetic spectra with 5 components (dashed), for different average energies (left) and widths (right panel). 
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5: Comparison between PL spectra and the cumulative 
sections for selected spectrum parameters. 
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FIG. 6: Reconstruction of the incident SN;/ spectrum from 
cross sections folded with a synthetic spectrum (dots). The 
dotted and full lines represent the 90% confidence levels for 
constructed spectra with 5 and 10% uncertainties on the pa- 
rameters a''''=i'^. 
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